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Abstract 

The squared eigenfunctions of the spectral problem associated to the 
Camassa-Holm equation represent a complete basis of functions, which 
helps to describe the Inverse Scattering Transform for the Camassa-Ifolm 
hierarchy as a Generalised Fourier transform. The main result of this work 
is the derivation of the completeness relation for the squared solutions of 
the Camassa-Holm spectral problem. We show that all the fundamental 
properties of the Camassa-Holm equation such as the integrals of motion, 
the description of the equations of the whole hierarchy and their Hamilto- 
nian structures can be naturally expressed making use of the completeness 
relation and the recursion operator, whose eigenfunctions are the squared 
solutions. 
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1 Introduction 



In this introductory section we shall give a brief account of the basic results 
about the spectral problem related to the Camassa-Holm equation (CH) [T^ . 
The CH equation 

Ut - Uxxt + '2ujUx + ZuUx - 2UxUxx - UUxxx = 0, (1) 

where w is a real constant, arises as a compatibility condition of two linear 
problems [12] 

^xx = (i + A(m + cj))* (2) 

= [j^-u^^x + Y'^ + l^ (3) 

where m = u — Uxx and 7 is an arbitrary constant. We will use the freedom 
provided by the presence of 7 for a proper normalization of the eigenfunctions. 

The CH equation models just like the Korteweg-de Vries (KdV) equation 
the propagation of two-dimensional shallow water waves over a flat bed. It also 
arises in the study of the propagation of axially symmetric waves in hyperelastic 
rods [32l [30] and its high-frequency limit models nematic liquid crystals [47l [9] ; 
more about the physical applications of this equation can be found e.g. in 
[I1[S1[53[31[351[^[11[Z5]. Both KdV and CH are integrable systems [551 U 
[Il[H[Tl[nT] (see also [H [31 HO]) , but while aU smooth data yield solutions of 
the KdV equation existing for all times, certain smooth initial data for CH lead 
to global solutions and others to breaking waves: the solution remains bounded 
but its slope becomes unbounded in finite time (see [HI [HI [1^ ) . The solitary 
waves of KdV are smooth solitons, while the solitary waves of CH, which are also 
solitons, are smooth if w > [3 IMl [701 [ZD [22] and peaked (called "peakons" 
and representing weak solutions) if w = [H [H [3 [3 [371 HI] ■ Both solitary 
wave forms for CH are stable [H UHl [SI] • The CH equation arises also as a 
geodesic equation on the diffeomorphism group (if lo — 0) [15l [23l [24l [60] and 
on the Bott-Virasoro group (if > 0) [571 [^ . 

The inverse scattering transform (1ST) for the CH equation and the related 
Riemann-Hilbert problem are considered e.g. in [T51 [201 [551 [55] . 

The 1ST can be realized as a generalized Fourier transform. The complete ba- 
sis of functions is represented by the squares of the fundamental ( Jost) solutions 
of the corresponding spectral problem. For the famous Zakharov-Shabat spec- 
tral problem (and its various generalizations) the problem is studied in detail and 
the results can be found in several important works, such as [Ml [501 I44j . 
The generalized Fourier transform for the Sturm-Liouville spectral problem, as- 
sociated to the fundamentally important KdV equation is also well studied and 
in this relation we can mention the books and articles [31 [HI [57l [36l HSl [58l [59] . 

The CH spectral problem ([2]) is of a weighted Sturm-Liouville-type. Our 
aim will be to construct the generalized Fourier transform which linearizes the 
CH equation as well as all equations of the whole CH hierarchy of integrable 
equations. Our main result is the derivation of the completeness relation for the 
squared solutions of ([2]) . We show that all the fundamental properties of the CH 
equation such as the integrals of motion, the description of the higher CH-type 
equations and their Hamiltonian structures can be naturally expressed making 
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use of the recursion operator L± and the completeness relation. In fact, the 
squared solutions are eigenfunctions of L± and the completeness of the squared 
solutions is the spectral decomposition of L±. 

In Section [2] are given all the necessary mathematical preliminaries about 
the spectral problem ([2]). This includes the definition of the Jost solutions and 
two sets of scattering data, as well as the time evolution of the scattering data. 

In Section [3] are presented the asymptotics of the Jost solutions for large 
values of the spectral parameter — ^ oo). The main difference with respect 
to the standard spectral problem, given by the Schrodinger equation, (e.g. like 
in the KdV case) lies in the fact that (l2|) is a weighted spectral problem, which 
requires different asymptotic expansions. 

The completeness relations of the eigenfunctions and the squared eigenfunc- 
tions of the spectral problem ^ are presented in Section [H The last one gives 
the possibility of expansion of an arbitrary function of the specified class over 
the complete basis. 

The Wronskian relations, derived in Section [5] allow to compute the gen- 
eralised Fourier coefficients for the potential of our spectral problem and its 
variation. 

The symplectic variables, i.e. the variables, expressed in terms of the squared 
solutions and satisfying the cannonical relations, with respect to a certain skew- 
product are given in Section [6l 

The recursion operator (whose eigenfunctions are the squared solutions) is 
computed in Section [71 

In Section [5] the whole Camassa-Holm hierarchy is constructed. The time- 
evolution of the scattering data for the hierarchy is obtained. 

In Section [9] the Hamiltonian structure of the Camassa-Holm hierarchy is 
explored. The canonical Hamiltonians are expressed both via the potential of the 
spectral problem and the scattering data. The hierarchy of Poisson structures 
and action-angle variables with respect to these structures is also obtained in 
this section. 

Finally, in Section (TU] the Inverse Scattering Transform for the CH hierarchy 
is outlined in the framework presented earlier for the CH equation alone |20j . 



2 Preliminaries 

In general, there exists an infinite sequence of conservation laws (multi-Hamiltonian 
structure) H„[m], n = 0, ±1, ±2, . . ., [HI [Ml Eal EH [39l [62] such that 

1 f 

Hi = — / mudx, (4) 

H. = l/(.» + ™» + 2.„=)d., (5) 

(a— a) — = (2Lua + mo + dm) — - — -. (6) 

dm dm 

The CH equation can be written as 

mt^{m,Hi}, (7) 
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where the Poisson bracket is defined as 

{A,B} = - f ^ {2Lud + md + dm)^dx, (8) 
J dm dm 

or in more obvious antisymmetric form. 

_ /(w + m)f— 9— —d—)dx (9) 
' J \Sm 5m 5m 5m J 

For simphcity we will consider the case where to is a Schwartz class function, 
w > and m{x,0)+Lij > 0. Then m{x,t)+LL: > for all t |16j. Let P = -j-Xuj, 
i.e. 

A(fc) = -i(fc^ + i). (10) 

The spectrum of the problem ^ under these conditions is described in [16] . 
The continuous spectrum in terms of k corresponds to fc - real. The discrete 
spectrum consists of finitely many points /c„ = n — I, . . . , N where k„ is 
real and < < 1/2. 

A basis in the space of solutions of ^ can be introduced by the analogs of 
the Jost solutions of eq. ([T]), f^{x, k) and f~^{x, k). For all real fc 7^ it is fixed 
by its asymptotic when a; 00 [12], see also [55] : 

lim e-*'=^/+(x,fc) = 1, (11) 

X — >oo 

Another basis can be introduced, f~(x,k) and/^(a;,fc) fixed by its asymptotic 
when X — > —00 for all real fc 7^ 0: 

lim e''=^/-(x,fc) = 1, (12) 

X — * — 00 

Since m{x) and oj are real one gets that if /^(x, fc) and f~{x, fc) are solutions 
of ^ then 

f+{x,k)^f+{x,-k), and r{x,k)=r{x,-k), (13) 

are also solutions of ([T]). The relations ([T5|) are known as involutions. 
In particular, for real fc 7^ we get: 

/±(a;,fc) =/±(x,-fc), (14) 
and the vectors of the two bases are related H: 

f'{x,k) ^ a{k)f+{x,~k) + b{k)f+{x,k), Im fc = 0. (15) 

The Wronskian W{fi, /2) = fidxf2 — f2dxfi of any pair of solutions of ^ does 
not depend on x. Therefore 

W^(/-(x, fc), /-(x, -fc)) = Wif+{x, -fc), /+(x, fc)) = 2ik (16) 

^According to the notations used in 1201 1211 f^(x,k) = ^{x,k), f~(x,k) = (p(x,k). 
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Computing the Wronskians W{f , /+) and V7(/+, / ) and using we 
obtain: 

a{k) = (2zfc)-iM/(r(a;,A:),/+(x,A:)), (17) 

b{k) = ~{2ik)-'Wif-{x,k)J+ix,-k)). (18) 

From and (fTB|) it follows that for real fc 



a(fc)a(-fc)-&(fc)&(-fc) = 1. (19) 

It is well known [16 that {x,k)er^^^ and f^{x,k)e^'^^ have analytic ex- 
tensions in the upper half of the complex fc-plane. Likewise f~^{x,k)e'^^^ and 
f~{x, k)e~'^^^ allow analytic extension in the lower half of the complex fc-plane. 
An important consequence of these properties is that a{k) also allows analytic 
extension in the upper half of the complex fc-plane and 

a{k) = a(-fc), h{k) = 6(-fc), (20) 

As a result (fTnj) acquires the form: 

Kfc)p-|6(fc)P = l. (21) 

In other words the relation (|17p is valid in the upper half plane Im fc > 
pniSO] , while (|18p makes sense only on the real line Im fc = 0. In analogy with 
the spectral problem for the KdV equation, one can introduce the quantities 
T[k) — a~^{k) (transmission coefficient) and TZ^{k) = b{±k)/a{k), (reflection 
coefficients - to the right with supcrindex (+) and to the left with superindex 
(— ) respectively). From ([21]) it follows that 



|T(fc)|2 + |7e±(fc)|2 = 1. (22) 



It is sufficient to know TZ (k) only on the half line k > 0, since from ([201 
a{k) = a(-fc), b{k) = b{-k) and thus n^{-k) = 1Z'^{k). Also, from 



Kfc)|2 = (l-|7^±(fc)n-^ (23) 

One can show that TZ^{k) uniquely determines a{k) fIT\ . 

At the points k„ of the discrete spectrum, a{k) has simple zeroes [16) . i.e.: 

a{k) = (fc - iK„)a„ + ^(fc - tK„)^a„ H , (24) 

and the Wronskian W{f^,f'^), pT|) vanishes. Thus and /+ are linearly 
dependent: 



/ (x.iKn) = bnf~^{x,iKn)- (25) 

In other words, the discrete spectrum is simple, there is only one (real) linearly 
independent eigenfunction, corresponding to each eigenvalue say 
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(26) 



From (f26|) and (fT2|) it follows that fni^) falls off exponentially for x 

±00, which allows one to show that fnix) is square integrable. Moreover, for 
compactly supported potentials m(x) (cf. ([251 and (fTKl) ! 

bn = b{iKn), 6(-lK„)=--^. (27) 

One can argue [68 , that the results from this case can be extended to Schwarz- 
class potentials by an appropriate limiting procedure. 
The asymptotic of /~ , according to ([Ti| , PT|) , P5)) is 



f-{x) = e^-^ + oie^-n, x^-oo; (28) 

f-{x) = &„e-''"" + o(e-''""), x^oo. (29) 

The sign of 6„ obviously depends on the number of the zeroes of f^. Suppose 
that < Ki < K2 < ■ • ■ < < 1/2. Then from the oscillation theorem for the 
Sturm-Liouvillc problem [S], /~ has exactly n — 1 zeroes. Therefore 

&„ = (-l)"->„|. (30) 

The sets 



5=^ = {7e=^(fc) (fc>0), K„, Ct = -^, n=l,...N} (31) 

tan 

are called scattering data. Throughout this work the dot stands for a derivative 
with respect to k and d„ = d(iK„), a„ = a(zK„), etc. In what follows we 
will show that each set - or S~ of scattering data uniquely determines the 
potential m(x). The derivation is similar to those for other integrable systems, 

e.g. [Miiiaiiiizs]- 

The time evolution of the scattering data can be easily obtained as follows. 
From ([TS]) with a; ^ oo one has 



lim (/-(x, k) - a{k)e-'^'' - h{k)e'^'') = 0. (32) 

X — *oo 

The substitution of k) into ^ with x ^ oo gives 

lim(/r-:^/-+7/") = 0. (33) 
From ([S^ . with the choice 7 = 7- ^ ik/2X we obtain 

at{k,t) = 0, (34) 
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Thus 



bt{k,t) 
a{k, t) = a(fc, 0), 



= |6(fc,i). 
h{k,t) = &(fc,0)e'^*; 



(35) 
(36) 



T{k,t) =T{k,0), 7^=^(fc,^) =7^±(fc,0)e=^^^ (37) 
Similarly, one can substitute f^{x, k) as x — oo 

lim (/+(x, k) - a{k)e'^^ + b{-k)e-'^'-') = (38) 

X — ^oo 

into Then the choice of the constant is 7 = 7+ = —ik/2\ and the final 
result (l36|) is, of course the same. 

In other words, a(fc) is independent on t and will serve as a generating 
function of the conservation laws. 

The time evolution of the data on the discrete spectrum is found as follows. 
iKn are zeroes of a(fc), which does not depend on t, and therefore Kn,t — 0. 
From ([27|) and ((35|l one can find 

C^,tW = C±(0)exp(±^i^t). (39) 

3 Asymptotics of the Jost solutions for |/c| ^ oo 

The analyticity properties of the Jost solutions and of a{k) play an important 
role in our considerations. We will need also the asymptotics of the Jost solutions 
for |fc| ^ 00 which have the form [201 HI] 



f^{x,k) = e^'=^+^'=/<i(\/'3fe)/"~i)''s' 



. nix) 



(40) 



where, for simplicity q{x) = m{x) + lj. 

An immediate consequence of the above formulae and (|17p is: 



lim a{k)e'°''' = 1, kc C+, 



where 



-Udx 



(41) 



(42) 



(43) 



Since a(fc) is independent, then a, as well as all the coefficients Ik in the 
asymptotic expansion 



00 J 

lna(/c) = -ia/c + ^-^. 



(44) 
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must be integrals of motion. The integral a is the unique Casimir function for 
the CH, |56| . One can easily check that {to, a} = 0, for the uniqueness argument 
and for the geometric interpretation see |56| . 

The densities, Ps of 1^ ~ J_^Psdx can be expressed in terms of m{x) using 
a set of recurrent relations obtained in [73l [511 [21] . 

Using the analyticity properties of a(fc) one can prove that it satisfies the 
following dispersion relation (fc G C+): 



liia(k) = -lak+y In- / ^ ^ ^ ^' ' dk' . (45) 

n—1 ^ 

where iKn are the zeroes of a{k). Thus we are able to recover the function a{k) 
in its whole domain of analyticity C+ knowing just its modulus |a(fc)| or the 
reflection coefficient TZ'^{k), for real fc, and the location of its zeroes, see 

The dispersion relation allows one to express the integrals of motion 
also in terms of the scattering data [21] : 

Is = - \n{l-\n^ikr)k'^dk-J2 ''\\^" . (46) 

These are known as the trace identities. In addition the integral a is expressed 
through the scattering data as follows. Note that for k — i/2, A(z/2) = from 
(fTOj) . In this case therefore the spectral problem ([2]) does not depend on to, and 
the eigenfunctions are equal to their asymptotics: f^{x,i/2) = e^5. From (fT7|) 
we obtain a{i/2) — 1 and from (|l5)l for k = i/2 we have 



a = li^ 

n=l 



/ 1 + 2^X^1 fMid^W),^. (47) 

Vl-2K„y TT Jo 4fc2 + l 



4 Completeness relations for the Jost solutions 

Here we outline the spectral properties of the linear problem 
Let us consider the function 



n/ M - / i^^k)f+{y,k) f {y,k)f+{x,k) 

Ri{x,y,k) = — e{y-x) + e{x - y) 48) 

ayk) ayk) 

where 9{x) is the step function. 

Lemma 1. i) Ri{x,y, k) is an analytic function of k for k G C+z 

ii) Ri{x, y, k) has simple poles for k = iKn', 

iii) Ri{x,y,k) is a kernel of a bounded integral operator for Imk > 0. For 

Imk — 0, Ri{x, y, k) is a kernel of a unbounded integral operator. 

iv) Ri{x,y, k) satisfies the equation: 
^ ' + Xqix)) Ri{x,y,k) ^4ikS{x~y). (49) 



dx^ V 4 
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Proof, i) and ii) are obvious. 

iii) For Imfc > the statement follows from the definitions pT|) and ([T] 
of the Jost solutions, which ensure that Ri{x, y, k) falls off exponentially for all 
x,y ±00. The same arguments for hnk = can only ensure that Ri{x, y, k) 
is a bounded function for x,y —> ±00. 

iv) is a consequence of ©, P7)l and the fact that 

d0{x — y)/dx — S{x — y). 

□ 

Remark 1. The function Ri{x,y,k)/{-iik) is a kernel of the resolvent for the 
linear problem 



The explicit expression for the resolvent PS]) can be used to prove the fol- 
lowing 

Proposition 1. The Jost solutions of 0) satisfy the completeness relation: 



-Ax -y) = ^r l^i^iM^dk + f2 c^f^ix)ftiy) (50) 

By dk is meant limi^oo /^^ dk and f^{x) = f^{x, inn). 
Proof. Consider the contour integral 

Ji{x,y)^7r-f Ri{x,y,k). (51) 

where the contour 7+ is shown on Fig. [TJ 

From the Cauchy residue theorem it follows that: 

N N 

:7i(a;,y) = ^ Res i?i(x, fc) = z ^ C±/±(x)/±(2/). (52) 

n—1 n—1 

In the evaluation of the residues we made use of (f2F]l and of the fact that 
e{x-y) + e{y-x)^l. 

Next we evaluate the integral J\{x^ y) by integrating along the contour. For 
the integration along the infinite semicircle we need the asymptotic of i?i {x, y, k) 
for |fc| 00. From and (jH]) we get 

Ri{x,y, k) - e^^-^^" ^""^ \(^) (^) ' + o(l/fc)l 9{y ~ x) + {x ^ y) (53) 



q{x) 



Only the leading terms in (j53p . which are entire functions of k contribute to 
the integral. This allows us to deform the infinite semicircle until it coincides 
with the real fc-axis. Then the integration over k is easily performed with the 
result: 

Ji.Ux.y)^ ^5{x-y). (54) 
q[x) 

To evaluate the integral over the real axis J^i^r{x, y) we will use the the fact 
that Ri{x, y, k) can be written in the form 

Ri{x,y,k) = lZi^^l^t^ + N,{x,y,k)9{x~y), (55) 
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Figure 1: The contours 7± = M U 7±.oc of integrations. 

where 

m{x, y, k) = ^[/-(x, k)f+{y, k) - }+{x, k)r{y, k)] (56) 

is an odd function of k and does not contribute to J^i^r{x, y). Indeed, from (jlSp 
we have 

N^{x,y,k) = L±l^il^.f+^^^k)[f+(y,-k)+n{k)f+{y,k)] 

f-{x,k)f+{y,k) + , Jf-ix,-k)-b{-k)f+ix,-k)] 
a{k) ^ a{~k) 

- TZ{k)f+{x,k)f+{y,k) 
rix,k)f+iy,k) f+{y,-k)f-{x,-k) 

a(k) a{—k) 

- TZ{-k)f+{x, ~k)f+{y, ~k) - TZ{k)f+{x, k)f+{y, k). 

Now it remains to equate the two expressions for 

Ji(a;,?;) = ^^ Ri{x,y,k)^Ji,oo{x,y)+Ji^R{x,y). (57) 
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to obtain the completeness relation for the Jost solutions. □ 

Remark 2. Another way to prove the completeness relations for the Jost solu- 
tions is to note that Ri{x, y, k) is a kernel of the resolvent defined for Imk < 
and then apply the contour integration method to the integral 

= dk Ri{x,y,k) - <j> dkRi{x,y,k). (58) 

5 Wronskian relations and generalized Fourier 
series expansion 

An important tool for the analysis of the mapping between the potential m(x) 
of eq. ([H) and the scattering data, related to this potential, are the so-called 
Wronskian relations. 

For what follows, we define the squared eigenfunctions 

F±(x,A:) = = zk„)- (59) 

Proposition 2. Let f{ki,x) and g{k2,x) he two eigenfunctions of the spectral 
problem @). Then the following identity (Wronskian relation) holds: 

I \ (I. \^ fx{ki,x)g{k2,x) - fiki,x)g^{k2,x) °° 

q[x)f(ki,x)g{k2,x)dx = —— (60) 

, A(fci)-A(fc2) 



Proof. It follows immediately from the fact that f{ki^x) and g{k2^x) satisfy 
Corrollary 1. From \6(J\) . {14^ , 177]) . J75)) . and from the formula 

lim P — ; — = TTi6{k), 

X — >oo k 

where P means Principal Value, one can obtain various 'orthogonality ' relations 
for the eigenfunctions, such as (cf. fMHi) 



q{x) f^ {ki , x) f {k2,x)dx = —2nuja{ki)5{ki — k2), (61) 
q{x)f^{x)f~{x)dx = iujdnSnp. (62) 

□ 

Let us define a scalar product as usual by 

/oo 
f{x)g{x)dx, (63) 
-oo 

provided the corresponding integral exists. 
We need also the skew-symmetric product 



11 



= / q{x){fx{x)g{x) - g^{x)f{x))dx, 

<J —OQ 

related to the Poisson bracket © via 



{A,B} = 



6A 6B_ 
5m ' Sm 



(64) 



(65) 



Some other useful Wronskian relation are formulated in the next two propo- 
sitions. 

Proposition 3. Let f{k, x) be an eigenf unctions of the spectral problem 
Then the following identity holds: 



{mx,f^) = ^\{fx - ffxx 



(66) 



Proof. Since f{x) satisfies ([2]), 

{fx f fxx) - 



{fl- ffxx)xAx = -\{mxJ^). (67) 



□ 



The right hand side of (|67p can be expressed also through the skew-symmetric 
product 



(m^, /^(x, k)) = {q^, /2(x, k)) = 



q{x) 



- 1 



(68) 



Corrollary 2. From Proposition[3[ and equations OI|J . i f5S|J and (SS) it follows 



{m,,F^{k)) 



F±(x,fc), 



q{x) 



-—a{k)b{Tk). (69) 



and t/ie limiting procedure gives on the discrete spectrum 



{m,,F^)^hq{x),F^ix)]]^0, (rUx, F^) ^ hqix), F^ix)]] ^ 

(70) 

where A„ = A(zk„). 

Another type of Wronskian relations relate the variations of the potential 
m(x) with the variation of the scattering data [21]: 



{f{x,k)bfx- fx5f{x,k))\ 



X — — 00 



dxXSq{x)f^{x,k), (71) 
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where Sf{x, k) is the variation of the Jost solution f(x, k) corresponding to 
the variation 5q{x) of the potential. The right hand side of ((7T|) can also be 
expressed through the skew symmetric product as follows: 



{f{x,k)5f^- fJf{x,k))\Z^ 



-\lf{x,k),5Q±{x)^, (72) 



where 5Q±[x) = [1/ \Jq[x) j^^Ay 5^/qijj). For the derivation of (17^ we per- 
formed an integration by parts and used the condition 5Q±{^oo) = ±5a = 0. 
From dni) and ^ we obtain: 



lF^{x,k),5Q±^ = T-fa\k)5n^{k), 



[[i^,?(x),<5Q±]] = 



^nCn 



aiSC- 



(73) 
(74) 

(75) 



Using ([7T|) one can also derive the following relations for the variations of the 
scattering data, for details see [U]: 



5a{k) 
Sm(x) 

6b{k) 
Sm(x) 
6 In A„ 
dni{x) 



2ik 



f+ix,k)f-{x,k) 



= —f+{x,-k)r{x,k) 

Ujbndn 



(76) 
(77) 
(78) 



The equations ^ - ^ lead further to 



sn^jk) 

5m{x) 

Sm(x) 

Sk„ 
5m{x) 

Sm{x) 



± 



m 

2ika'^{k) 

An 



F^{x,k) 



1 



\4:iuJKnXn 



^\f^{x)+F^{x) 

Our J 



-KCtFt^ix) 



(79) 
(80) 
(81) 
(82) 



These relations allow us to calculate the skew symmetric products between 
the squared solutions. 

Proposition 4. Let fi{ki,x), f2{k2,x), gi{ki,x), g2{k2,x) be eigenf unctions 
of the spectral problem (0). Then 



[[/l(fcl)5l(fcl),/2(fc2)ff2(fc2)]] = 



{.fidx.f2 - f2dxfi){9idxg2 - g2dxgi)) 



A(fci)-A(fc2) 



X — — OQ 



(83) 
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Proof. Straightforward to check using the fact that all functions are eigenfunc- 
tions of the spectral problem - see [H] for details if necessary. □ 

Corrollary 3. From Proposition^ and hS^) . \38\) it follows 

f/T(fc2)]] = ±27rwA(fci)5(fci - k2), kia G K (84) 
where U^{x,k) = 2ik!i^k) ^^(^^ diserete speetrum 

, i^,^ ]] ^ ^2»c.«„a2 ^„„, , [[F+ , F~]] = 0, , = 0. (85) 

Up to now we demonstrated that the mapping between the scattering data 
TZ^{k) and the potential m{x) (or, equivalently, q{x)) is expressed through the 
squared solutions of ([2|). The same squared solutions relate also the variations 
of the reflection coefficients STZ^{k) with the corresponding variations of the 
potential 5m{x) (or, equivalently, 5q{x)). 

In order to ensure that these mappings are one-to-one we have to prove that 
the squared solutions form a complete set of functions in the space of allowed 
potentials; in other words we will prove that the functions F^{x, k), F^{x) form 
a basis in the space of allowed potentials. To this end we consider the function 



F-{x^)F+{y^ 

R[x,y,k)^ ^^2(fc-) Gyy-x) (86) 

, 2/+(x, k)f-{x, k)f+{y, k)f-{y, k) - F-{y, k)F+{x, k) _ 
+ ka^{k) ^^"^ y> 

Lemma 2. i) R(x,y,k) is an analytic function of k for k € C+; 

ii) R{x, y, k) has second order poles at k — iKn; 

iii) R{x, y, k) is a kernel of bounded integral operator for Irak > 0. For Imk = 

Ri{x,y, k) is a kernel of a unbounded integral operator. 

Proof, i) and ii) are obvious. 

iii) For Imfc > the statement follows from the definitions pT|) and 
of the Jost solutions, which ensure that R{x,y, k) falls off exponentially for all 
x,y ±oo. The same arguments for Imfc — can only ensure that R{x, y, k) 
is bounded function for x, y — *■ ±cxd. □ 

Proposition 5. The following completeness relation holds: 
-.eix - y) = / dk 



^q{x)q{y) 27ri ka'^{k) 



Y.-—^\^n{^)Fi{y) + F-Ax)F+{y)-(— + '^)F-{x^^^ (87) 

n—1 "' 
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Proof. The proof is similar to the one of Proposition [TJ 
Consider the contour integral 

Jix,y) = TT^ / Rix,y,k)dk. (88) 

where the contour 7+ is shown on Fig. [TJ From the Cauchy residue theorem it 



(89) 



follows that: 






N 

- Res R{x, y, k) 

n=l 


N 

= E 

n=l 


-l^\F~{x)F+{y) 



^~{x)F+{y) + F-{x)F+{y) -(— + -)F-{x)F,t{y) 



In this evaluation we used (^5)) and the fact that 9{x — y)+9{y — x) — 1. 

Next we evaluate the integral J{x,y) by integrating along the contour. For 
the integration along the infinite semicircle we need the asymptotic of R{x, y, k) 
for |fc| 00. Due to ([40|) and ([IT]) we get 

R{x,y,k) = I e{x-y){l + o{l/k)) (90) 

+ i (Roix, y, k)e{y - x) - R^{x, y, -k)e{x - y)) (1 + o{\/k)), 



Ro{x,y,k) = " 2..gd.V^. 

Only the leading terms in (|53p , which are are entire functions of k contribute 
to the integral. This allows us to deform the infinite semicircle until it coincides 
with the real fc-axis. Then the integration over k is easily performed with the 
result: 

J^ix, y) = —^==e{x ~ y). (91) 

To evaluate the integral over the real axis J^R{x,y) we will use the the fact 
that R(x, y, k) can be written in the form 

= nf.iElMl.^!li£^e,.-y) ,92) 

where Ni{x, y, k) is defined in ((55)) . The second term in ([55]) is an odd function 
of k and does not contribute to jR(x,y). 

Now it remains to equate the two expressions for 

^^^^y) = ^ij^ Rix,y,k)dk = Jo,{x,y) + jR{x,y). (93) 

to obtain the completeness relation for the squared solutions. □ 

Corrollary 4. The completeness relation can be rewritten in the following 
equivalent form: 
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JV 1 ■■ 

+ T.-—^\^n{^)F:i{y) + F-{x)F+{v) - (— + ^)F~{x)F^{y) . 

where 

F^{x, k) = J ^F^{x, k). (95) 

V LO 

The completeness relation allows one to expand any function X{x) over the 
squared solutions F^{x, k) or F^{x, k). To this end we multiply both sides of 
([57)1 with ^qyX{y) + q{y)Xy (resp. with ^qxX{x) +q{x)Xx) and integrate over 
Ay (resp. dec). A simple calculation gives: 



(96) 

provided X{x) vanishes for x —> zLoo. The expansion coefficients ^^(fc) are 
given by: 

e±(fc)= dkF^{y,k) Qg,X(y) + g(y)X(y)^ -i[[F±(a;,fc),X(x)]], 

(97) 

where [[•,■]] is the skew-symmetric product . Similarly for „ and „ 
we find: 

^lu = ^ ^(^)]] , = I (^), x{x)]] - (-^ + 1^) ci,„, 

(98) 

Analogously, if we multiply the completeness relation by ATj, (resp. by 
Xx) and integrate over dy (resp. dx) we get: 



N 

- f PT^± j_ pTf 

X,n) ' 

(99) 

where the expansion coefficients C^(fc)j ^-nd are expressed by another 
skew-symmetric product: 

/OO 1 
dfci?±(2/,fc)X,(y) = -[[F±(a;,fc),X(x)]], (100) 
-00 A ~ 

Obviously the two skew-symmetric products are related by a gauge-like 
transformation. Indeed, if F^{x,k) and F^{x,k) are related by ([5^ and if 
X(x) = \/(7(x)/u;X(x) then: 

lF^{x,k),X{x)]]^iu[[F^{x,k),X{x)]]. (101) 

The following lemma demonstrates that the mapping between X{x) and its 
set of expansion coefficients is one-to-one. 
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Lemma 3. i) A necessary and sufficient condition for X{x) to vanish is that 
all its expansion coefficients vanish: 



ef(fc)=0, ei„ = 0, ef,„=0. (102) 
ii) A necessary and sufficient condition for X{x) to vanish is: 

el(fc) = 0, ||„ = 0, e|„ = 0. (103) 

Proof, i) Let us assume that X{x) — for all x. Substituting it into the right 
hand side of ^ we immediately get (fTU^ . Let us now assume that ()102p 
holds. Inserting it into the right hand side of the expansion wc find that 
X{x) = for all x. ii) is proved analogously using (jlOOp and the expansion 

(uni). □ 

Remark 3. We introduced two sets of squared solutions F^{x, k) and F^{x, k) 
related by the gauge-like transformation Since q[x) is always positive the 

two sets are obviously equivalent. 

6 Symplectic basis and canonical variables 

Here we introduce a third set of squared solutions, known as symplectic ba- 
sis whose special properties will become clear below. The symplectic basis was 
introduced for the first time in analyzing the nonlinear Schrodinger (NLS) equa- 
tion in [43j. 

The canonical (action-angle) variables for the CH equation are known pT] : 

9k k 
P^^^ ^ ^:7i:\ll^l«(^)l= ^7TTJln(l-7^±(fc)7^±(-fc)), (104) 

1 V + lk) 

m EE arg6(fc) = -ln^^, (105) 

Pn = A-\ 0„=ln|6„|, (106) 
satisfying 

{p(fci)Xfc2)} = 5{ki-k2) (107) 

{p{ki),p{k2)} - {0(fci)Xfc2)} = O, A:i,2>0 (108) 

{Pm,4>n} Smn, (109) 

{Pm,P«} = = 0. (110) 

The following proposition gives the symplectic basis in the space of Schwartz- 
class functions: 

Proposition 6. The quantities 

ZmLoMk) \ / 



J— (n+ik)F+ix, k) - n+{-k)F+{x, -k)) , (111) 



2'KiLuX{k) 
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Q{^,k) = -—^^^^(n-{k)F-{x,k)+n+{k)F+{x,k)), (112) 

1 . 1 



Pnix) = — C+F+(.t) = — C-F-(a;), (113) 

Qn(:r) = ^{c+F+{x)-C-F-ix)). (114) 

satisfy the following canonical relations: 

[[P{ki),Q{k2)]] = <5(fci-fc2), (115) 

lv{ki),V{k2)]] = [[Q(fci), Q(fc2)]] = 0, fci^2>0 (116) 

[[/'m.On]] = S„,n, (117) 

[[-Pm,P„]] - lQm,Qn]]=0. (118) 

Proo/. Due to §5^, 



om[x) dm[x) 

and with ((711) - ([Hll) we obtain pTTjl - (fTT^ . From (fTUS)) . (fTUH)) one can see that 
(j)n can be obtained by a Umiting procedure fc iHn from 0(fc) (up to an overall 
multiplier i). Therefore, Qn{x) can be obtained in the same way from Q{x, k), 
giving (dm). Note that PT7)) is satisfied due to pT^ and (fTTI)l . □ 

The usefulness of the symplectic basis (jllip - (|113p is more evident from 
the following two propositions. 

Proposition 7. 

—l===s{x-y) = / (v{x,k)Q{y,k)~Q{x,k),V{y,k))dk 

+ ^ (p„(a:)Q„(y)-g„(x)P„(2/)). (121) 



where eix — y) = — y) — 6{y — x)^ /2. 
Proof. From (jllip we have the identity 

'R-{k)F-{x,k)+n+{k)F+{x,k) = 

n-{-k)F-{x,-k) + n+{-k)F+{x,~k). (122) 

Applying several times p22p one can check further that 



Aoj(p{x,k)Q{y,k) - Q{x,k)P{y,k) 

1 ^ F-{x,k)F+{y,k)-F+{x,k)F-{y,k) 
2m\ ka^{k) 



+ (fc->-fc)) (123) 
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With pi3p and (|114p it is straightforward to obtain 



N 

4cjY, {Pn{x)Qn{y) - Qn{x)Pn{y) 



f(^"(--)^»"(^)-.y(-)^"(^)+(x <->,)). (124) 



n—1 

Now, (fT^ foUows from and the representation dST]). □ 

The completeness relation (I119p allows one to expand any smooth function 
X[x) vanishing for x — > ±00 over the symplectic basis. To this end we multiply 
both sides of pi9p by ^qyX{y) + q{y)Xy and integrate over dy. The result is: 

X{x)^ / dk{V{x,k)dpx{k)-Q{x,k)px{k)) + Y,iPn{x)(Pn,X-Qn{x)pn,x), 

(125) 

where 



<t>x{k) - lQ{v,k),X{y% px{k)=p>{v,k),X{ 
K,x = \Qn{y),X{y% p„,x= [[P«(?/),X(2/)]]. (126) 

Corrollary 5. The completeness relation 1121]) can be cast in the following 
equivalent form: 



^six-y)= {v{x,k)Q{y,k) ~ Q{x,k),V{y,k)yk 

N 

[Pn{x)Qn{y) - Qn{x)Pn{y)) , (127) 



where the elements of the new symplectic basis V{x,k), Q{x,k) are related to 
the old ones by: 



f{x, k) = Vq{x)r{x, k), Q{x, k) = y/q{x)Qix, k), (128) 
for all k eM.(J {iHn}. 

With (|127p the analogue of the expansion (|125p for any smooth function 
X{x) vanishing for x — > ±00 is: 

l(x)= / dk(v{x,k)4>x{k)-Qix,k)px{k)) 

N 

+ ^ (A^^X,n - Q«PX,„) , (129) 

where the expansion coefficients px{k), px,n and (f>x{k), (f>x,n are given by: 

px{k) = [[Q{k),X]], 4>x{k)^[[V{k),X]\, (130) 



19 



~PX,n = [[Qn,^]], ^x(fc)= [[P„,l]], 

where [[ •, • ]] is the new skew-symmetric product (jlOip . 

The following lemma demonstrates that the mapping between X{x) and its 
set of expansion coefficients is one-to-one. 

Lemma 4. A necessary and sufficient condition for X{x) to vanish identically 
is that all its expansion coefficients vanish: 

Pxik)^0, 4>xik), Px,„ = 0, 4>x,n = 0- (131) 
The proof is analogous to the one of Lemma [31 

7 Recursion operator for the Camassa-Holm hi- 
erarchy 

As in Section [5] we can view the completeness relations (|87l) and (jM]) as spectral 
decompositions for the recursion operators L± and L± defined below. 

Proposition 8. Let us define the recursion operators L± and L± and their 
inverse L± and L± as follows: 

L±=V-^M{x)-2 dym'iy) . L±^—V~^ dy ^/^d,- , {132) 
L±^-^ f dy -^Vydy •, L± - / dy -^dyVy^ ■ .(133) 
where = 9^ — 1. Then the following relations hold: 



L±F^{x,k) 


= \F^{x,k), 


L±Ft{x) 


L±Ft{x) 


= ^Ft{x) + 


^'':;f.^{x). 

ujXi 


L±F^{x,k) 


- \F^{x,k), 


L±Ft{x) 


L±Ft{x) 


= ^F,t{x) + 




L±F^{x,k) 


= XF^{x,k), 


L±F^{x) 


L±Ft{x) 


= XnF.^ix)~ 


LjJ 


L±F^{x,k) 


^ XF^{x,k), 


L±Ft{x) 


L±Ft{x) 


= Xj^ix)- 





1 



(134) 



^ Ft{x), (135) 



Xn 



OJ 



(136) 
(137) 



Proof. It is not difficult to prove that L±L± and L±L± act as an identity 
operator on any function X{x). In order to prove (|134p one can make use of the 
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fact that it can be reformulated as follows: 



dl - \\qix) + 2A / dy™'(2/) . = 

±OC> 



(138) 



This follows from the fact that f^{x, k) are eigenfunctions of the spectral prob- 
lem ([2]) with the asymptotics (jlip and (fT2|) . The rest of the relations are proved 
analogously. □ 

Corrollary 6. The recursion operators L±, L± and their inverse L±, L± sat- 
isfy the following relations: 

[[L+X, y]] - [[X, L^Y]] , IL+X, r]] - [[X, L^Y]] , (139) 
[[L+X,Yl = [[X,L^Y]], [[L+X,Yl^[[X,L^Y]], (140) 



for any pair of functions X{x) and Y{x). 

Proof. It follows easily from the definitions of the skew symmetric product and 
of the recursion operators using integration by parts. 

□ 

Let us now define the operator L = ^(i+ + L-) cf. (|132|) . and L = ■^(i+ + 
L_). A simple direct computation shows that a kernel of L for w ^ is empty, 
therefore it is possible to define the inverse operator L~^. It is clear from piip . 
(Uni), (USUI that 

LVix, k) = \-^V{x, k), LPn{x) = X-^Pnix). (141) 

One can compute the action of L to the remaining part of the symplectic 
basis, making use of the following proposition: 

Proposition 9. The following relations hold: 

L^F^ix, k) = \f^{x, k) ± ^a(fc)6(Tfc). (142) 
A A 

Proof. One can write ()142p as 



dl-l- 4Aq(x) + 2A / dy m'{y) . {f^f T —a{k)b{Tk) = 0, 

which is fulfilled due to the fact that are eigenfunctions of the spectral 

□ 



problem (O with asymptotics ((TTl), ([Ti| . 
Corrollary 7. From the above Proposition and illS^) it follows 



L±Q{x,k) = -Q{x,k)±2k. 
A 



(143) 
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Now it is clear that the eigenfunctions for the operators L and L are the 

elements of the symplectic bases V{x,k), Q{x,k), Pn{x), Qn{x) and V{x,k), 
Q{x,k), Pn{x), Qn{x), e.g. 

LV{x, k) = jV{x, k), LQ{x, k) = jQ{x, k), (144) 

LPnix) = -^Pnix), LQn{x) = -^Qnix). (145) 

Proposition 10. The inverse of the recursion operators L and L are given by: 

L-' = L=^{L+ + L_), L-i = L = i(l+ + L_). 

Proof. Checked by direct calculation. □ 

Corrollary 8. The recursion operators L, L and their inverse L, L are 'self- 
adjoint' with respect to the skew symmetric product: 

ILX, y]] = [[X, LY]] , ILX, r]] = [[X, LYl , (146) 

\[LX, ¥]] = \[X, LY ]] , ILX, Y]] = \[X, LY ]] , (147) 

for any functions X{x) and Y{x). 

8 Expansions over the squared solutions and the 
CH hierarchy 

It has been demonstrated that the squared solutions satisfy the completeness 
relation and therefore can be considered as generalized exponents. In this section 
we will derive the expansions of three important functions and demonstrate how 
they can be used for establishing the fundamental properties of the Camassa- 
Holm hierarchy. 

The first of these functions is the 'potential' m{x), or, rather one of the 

functions yjuj/q{x) — 1 or q{x)/LU — 1 , which are completely determined by 
m{x) and vice-versa. Its generalized Fourier coefficients are determined by the 
scattering data. 

Proposition 11. 

-2 / kV{x,k)dk + 2y2i^nPn{x), (150) 

•^0 n=l 
/.oo N 

-2 / kP{x,k)dk + 2y2i^nPn{x). (151) 

•^0 n=l 



q{x) 
1- 



q{x) 
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Proof. The first expansion (|148p is obtained by multiplying both sides of ([57)1 
by ruy, integrating with respect to y and using ((^ . ((701) . The second one p49p 
follows from the first one and from (^5)) . The expansion coefheients of the third 
expansion (|150p can be calculated using ()126|) . the definition of the symplectic 
basis (|lll|) - (|114p and ((Ml), COl). The fourth expansion p5ip is an immediate 
consequence of the third pSOp and (|128p . □ 

The generalized Fourier expansion for the variation of the 'potential' reads 
as follows: 

Proposition 12. 

5Vq{y)dy = — / -—5Tl^{k)F^{x,k)dk (152) 



^/q(x) J±oo 27ri \{k) 



^ r 1 r± 

± E IT - ^nS\n)Ft (X) + -^S^nF^ (x) 
LAn I An 

n—1 



±oc> 



±CXD 



SJ^dy = ^ r -!—5n^{k)F^{x,k)dk (153) 
V 27^^ A(A:) 

^ r 1 r± • n 

± E IT ^^^^ - ^n6\n)Ft (X) + -^S^nF^ (x) 

^ LA„ iA„ J 

n— 1 

5\l^^dy = dk{v{x,k)5<i){k)-Q.{x,k)5p{k)^ 

N 

+ ^(P„(a;)<50„-Q«(a;)^p„). (154) 

n=l 

Proof. The expansion (|152p follows from (|96p choosing v'C(a;) = 5Q±{x). The 
corresponding expansion coefhcients (|97p. ([M]) are expressed in terms of the 
scattering data variations using ([751) -([75 1) . Eq. ()153p follows from ()152p and 

([55)) . The expansion (fTM| follows from (fT^ with l(a;) = /|_^ S^J^dy. Note 
that condition = ensures that the left hand side of (|154p is independent of 
the choice of the lower limit of the integration. The corresponding expansion 
coefficients are evaluated using the definition of the symplectic basis (|llip - ()114p . 
the Wronskian relations (I731)-(17S1) and 

□ 

The expansions p52p - p54p are valid for all variations of the potential 5m{x) 
preserving the value of the integral a. An important subclass of these variations 
are due to the evolution of r7i(a;, i). 

Effectively we consider a one-parameter family of spectral problems, allowing 
a dependence on the additional parameter t, such that m{x, t) is a Schwartz class 
function for all values of t. The variation of the potential with respect to t is 
given by: 

5m{x, t) = m{x, t + St) - m{x, t) ~ mtSt + 0{{5tY). (155) 

For such potentials the corresponding scattering data, e.g. TZ'^^k.t), C^{t), Kn 
generically will depend also on t. Keeping only the first order terms with respect 
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to dt we find that the corresponding variations of the scattering data are given 
by: 

sn^{k) = nt{k)st + oiistf), sc^ = c^^^st + oastf). (i56) 

With all these explanations from Proposition [T^ one easily proves the following 
Corrollary 9. 

/ iV^)tdy = — / -^nf{k)F^{x,k)dk (157) 
Vq{x) J±oo 27ri A(fc) 

LA^ I An -I 

71—1 

71—1 ^ 

dy = dk{v{x,k)Mk)-Q{x,k)pt{k)^ 

N 

+ ^ (P„(x)0„,t-Q„(x)p„,t) . (159) 



±oo 



±oo V ' / i -^0 



Proposition 13. Let ri(z) be a rational function such that its poles lie outside 

N 

the spectrum U Then: 

n—l 

n=l 



g(a;)N, _ 1 /""^ 2fc7^±(fc) ^ ~± 



n{L) 



+ E^^n"(^"')^n(^)' (161) 
n—l ^ 



1 ) = -2 / kn{\-^)V{x, k)dk + 2 V K„17(A-i)P„(x),(162) 



g(a;) 



f](L) n - j ="2jJ^ m(A"i)P(a;,fc)dfc + 2^K„l](A~i)P„(a;).(163) 
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Proof. Eq. (|160p is obtained by acting on the expansion (|148p with the operator 
ri(L±) and using (|134p . The condition imposed on the poles of n(z) ensures 
that ri(A~^) and r2(A,7^) are all finite, so the right hand side of (|160p is well 
defined. Eq. (fTCT|) follows from (fTig]) and (fHS]) . Eqs. (flH^ and (fTH51) are 
derived analogously using the expansions (|150p and ()151|) and ()144|) . 

□ 

Corrollary 10. //m Proposition \13\we use the operators n{L±), il{L±), n{L) 

and n{L) respectively, then in the right hand sides of il60\) - ri63\) the factors 
r2(A~^) and rt(X^^) will be replaced by il{X) and ri(A„) respectively. 

Proof. Follows easily from the arguments used in Proposition 1131 and (|136p and 
(IWl) . □ 

Now it is easy to describe the hierarchy of Camassa-Holm equations. To 
every choice of the function ^{z), known also as the dispersion law we can put 
into correspondence the Nonlinear Evolution Equation (NLEE): 

^£(V^).d, + f7(£)(y|-l)^0, (164) 

£ iV^)tdy + n{C) (l - ^i) = 0, (165) 

where £ (resp. C) is any of the operators L_ or L (resp. Z+, Z/_ or L). 

What will be demonstrated below is that the hierarchy (|164p (resp. (|165p l 
can be generated by each of the recursion operators L±, L, (resp. L±, L) 
and their inverse. The fact that the expansions over the squared solutions 
provide the spectral decompositions of the recursion operators makes evident 
the interpretation of the ISM as a generalized Fourier transform. Using these 
expansions we will show that each CH-type equation is equivalent to a linear 
evolution equation for the scattering data. 

Proposition 14. i) Each of the NLEE {164^ and 11165]) is equivalent to the 
following linear evolution equations for the scattering data: 

7^f(fc)T^fc^^(A"^)7^=^(fc) = 0, (166) 
Ct,t±^n^{K')Ct ^0, (167) 

Kn,t = 0. (168) 

ii) Each of the NLEE {164^ and hl65\) is equivalent to the following linear 
evolution of the action-angle variables: 

Mk,t)-kn{\-^) ^0, pt(fc)=0, (169) 

+ K„f^(A-l) = 0, Kn,t = 0. (170) 

Proof, i) Let us consider the NLEE p64p fixing up £ = L± and let us expand the 
left hand side of (|164p over the squared solutions {x, k, t) using the expansions 
(I157P and p60p . It is easy to check that the corresponding expansion coefficients 
coincide with the left hand sides of (1166p - ()168p . It remains to make use of i) in 
Lemma [3l To prove the equivalence of (|165p to the linear eqs. ()166p - (jl68p we 
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fix up £ = L± and make use of the expansions (jl58p . (|16ip and of ii) in Lemma 

131 

ii) Now we choose C = L and expand the left hand side of (|165p over the 
symplectic basis using ()169p . (|170p . It remains to apply Lemma HI □ 

The NLEE (fTBi]) and ([1^ can be simplified as follows: 

'?t + ^/g[^/g^^(/:)(y|-l)]^ = o, (iri) 

,.+ [l7(£)(l-y5)]^ = 0, (172) 
Example 1. FFif/i ri(^) — z one can easily check that 

L±(^^-l)^2u (173) 

and thus the equation ([177]) becomes the C amass a- Holm equation ([I]). 

The higher degree polynomials in (|17ip produce the other members of the 
Camassa-Holm hierarchy. All members of the hierarchy share the same spectral 
problem ([2]), and thus their solutions have the same x-dependence. The only 
difference is the time-evolution of the scattering data of the members of the 
hierarchy. 

Now it is clear how to extend the dispersion law il{z) for the case of ratio of 
two polynomials: 



Corrollary 11. Let 



where fii(z) and ^2{z) are two polynomials. The corresponding NLEE can be 
written in the form: 

ni{L±)^£ (V9)t% + f^2(i±)(y|-l) =0 (175) 

is also equivalent to \16b]) - hl68\) where il(z) is given by {174^ . 

Example 2. ^1(2;) = z and Vl2{z) — 1, i.e. Vl{z) — Xjz. The equation |J75[ ) 
due to USS^) has the form 



qt+dM-'^)^'^^0, (176) 

which is exactly the extended Dym equation Jl^. \40[ \ 69^ . 

Of course, further generalizations are possible e.g. by introducing another 
time-like variable, see '11'. 
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Corrollary 12. The functional derivative can he expanded over the sym- 

plectic basis as follows: 

Sf 



5m{x) 







dk 



({f,m}V{x,k)-{f,p{k)}Q(x,k)] 

N 

+ E ({/' 'I'r^Pnix) ~ {/, Pn}Qn{x)) . (177) 



Proof. Insert X{x) = into cq. (fT25l) and (fTM]) . In view of (fTTOl) and 

(|120p and (p5)) all expansion coefficients turn into the corresponding Poisson 
brackets. □ 

Corrollary 13. The Poisson bracket {/, 5} can be expressed as: 
{/, 5} = ({/, Pik)}{9, m] ~ {/, m}{9, Pik)}) dk 

N 

+ E ({^' P"^^9, M - {/, Mig, Pn}) . (178) 

n=l 

Proof. Multiply both sides of (|178p by j^^^ and take the skew-symmetric prod- 
uct. Due to 1^5]) the left hand side becomes {/, 5}. The right hand side follows 
from (inil) and (HH]), (fHUl) . 

□ 

Corrollary 14. From [TT^ with f ^ Hi ^ ^ J mudx and the relations {0, Hi} = 
k/X{k), {(j)n,Hi} = —Kn/Xn, scc ^21] wc obtain 

<^)=- Y7VT'P{x,k)dk + J2^Pr,{x). (179) 

Jo ^i 

Since the potential m{x) has the meaning of 'momentum' [S71[5n], (|179p gives 
expansion of w = (1 — dx)~^m over the 'modes' 'P{x, k), P„(x), representing the 
'momentum' part of the symplectic basis. Actually, (|179p is equivalent to (|148p 
due to p73p and the fact that 

L^V{x,k) = \-\k)V{x,k), L±Pn{x) ^ X-^Pn{x), (180) 

see dm]), dm, (UHl). 

9 Hamiltonian formulation for the CH hierarhy 

Let us start from the following observation. From the identities (|173p . (|179p we 
have 



(181) 



2 vy 9 
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(182) 

and therefore, taking into consideration pSOp . (|14ip . 

«M/f-')-«w(/f-')- 

Thus, the equation (|171|) can be written as 

[[<5(.),il7(i)(y|~l)]], (183) 

where <5(») denotes the delta function. 

Due to (|65p we can write this equation in Hamiltonian form 



mt^{m,H^}, (184) 
with Hamiltonian such that 

From (|162p we have further 

1 j — - pOO ^ 

-n{L)( --i) =- / kn{x-^)v{x,k)dk + y2KM{K^)Pn{x) 

Jo om{x) 



N . 

1^ (186) 



5m{x) 

Therefore, in view of (fTSS]) and (fTSe]) 



m(A-i), = «„r!(A-i), (187) 



5p{k) ' dp 

and finally 



/■ 2 f K ^ 

h''^- kn{\-^)p{k)dk--y -^n{x-^)dK„. (188) 

Example: As an example we can point out the CH equation (il(z) = z). 
The expression (|188p gives 



„a_„ ri-Ui.^ 2V-/, l-2^„ , 4^„(l + 44) 

The last expression was obtained in a different way in pij.D 
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Noticing that = 9^. ^ we can rewrite ^ as 



dm 2 dm 



due to p32p . or as 



^ = --L — " ^ 190 
orn 2 om 

If H^^ = is the Hamiltonian with respect to the Poisson bracket ([S]), the 
other conservation laws can be generated according to 

dH^lm] ^ l^ dH^l,[m] ^ ^^^^^ 



5m 2 Sm 

For example 



6m 2 Sm 

N 



^ ^-^0 n=l 



or in general 



Sm 



n 







■^«„(-2A„)'"^f^(A-i)P„(a;); 



' - fc(-2A)i-^"f7(A-i), ^ = K„(_2A„)i-^'r!(A-i), (192) 



for J > 1, or 



AT 

Hf = - I k{~2X)'-m{X-')p{k)dk - - V / ^(-2A„)i--'r!(A-i)dK„ 

1193) 

Example: For the CH equation {fl{z) = z) the expression p93p giveo 



■^This is also the quantity given in formula (81) of [21) . however with a technical error in 
the contribution from the continuous spectrum. The correct expression should be 

3^/ 1-2k„ 4k„(3 + 32k2 -48k*)\ 2^uj^ [°° k'^ In \a{k)\ 



^^=- i^y-TT^^ — 3(i-4«^)3 — i+— /, 



(4fc2 + !)■ 



'-Ak. 
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^ n—1 



^ 1 - 2k„ 4k„ (3 + 32k2 - 48k^ 



3(1-4^2)3 

(194) 

□ 

Since L^^ is well defined, it is possible to consider the formal expression 
(fTMl for j < 0. 

Example: For the CH equation (|193p gives (j = 0, ri(z) = z): 



(195) 



In terms of q{x) we have0 

/oo 
(y^-^/ZJ)2da:. (196) 
-oo 

Indeed, since SHq / 6m — 1 — yjuj/q, 6 Hi /6m — u, Hq and Hi are related through 
(fT89ll . see ([T73| .a 

Example: With (|189p . (|132p one can check that for the CH equation 

On the other hand, for j — —1 and Q{z) = z, (|193p gives 

H.i^-A A:Ap(fc)dfc + 2 V ( In — ^ ^-4k„). (198) 

Jo _-, ^ 1 ^ 2k„ / 



□ 



CorroUary 15. Obviously, ^^(^)~^/^ — ^h!?^^^"^, therefore the expressions 
on the right hand side of j jff7| ), {iff<$|) give i/ie Hamiltonian of ^176^ , (up to a 
constant factor) with respect to the Poisson bracket 

Here we notice that for the CH equation {^{z) — z) the integral a, 
(071) 

°-£(y!-')^^^|:>"(l^)"+r^»)* <™' 

is not of the form (jl93p . However, the hierarchy is generated by a{k) and 
a = limfe^oo ^l^^'i(^)- Thus, a from (|199p does not give rise to a separate 
hierarchy {«„} of conservation laws. Another way of seeing this is to assume 
the contrary, that there exists an independent hierarchy, such that 



P^-L-^, n = ±l,±2,.... 

6m 6m 



*Note the difference in comparison with the definitions of the integrals in |21| 
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However, one can check that and 4^ are related as follows: 

' dm dm 

Sm 6m Sm ' 

and therefore no other independent integrals arise. 

Example: The following integrals are related to a and an integral from 
(fTMD : 



mdx = Hn + 2uja 



rip(fc)d. + 2.Vflni±^ + ^ 



/ L 



I-i 



etc. □ 

The equation ()190p can be written as 



g2 N 1 
-|-jdx = i/_i--a, 



Therefore, defining a Poisson bracket 

{Ai3}(n) 

we can write the equation (jl7ip as 



5ra 



— ,(-2L-i)"-i — 



(200) 



(201) 



(")■ 



(202) 



The following statement gives the canonical variables for (|171|) with respect to 



Proposition 15. Let us define 



p{k)^^) = p{k), 0(fc)(„) = (-2A)i-'V(A;) 



Pu 



(203) 
(204) 



These variables satisfy the following canonical relations with respect to the bracket 



{p{kl)(n)A{k2)(n)}(n) 
{p{kl){n),P{k2)(n)}(n) 
{Pl,(n),4'p,(n)}(n) 
{Pl,{n)-, Pp,{n)}(n) 



5{ki - k2), 

{'/>(fcl)(n),0(fc2)(n)}(n) = 0, fci,2 > 

(n), 0p, («)}(«) = 0. 
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Proof. One can perform the verification making use of p44p . (|145p . For the 
quantities on the continuous spectrum the computations are particularly simple, 



e.g. 



{P{kl)(n),(f>{k2)(n)}{n) = (-2A(fc2))^ 



-2A(fc2))i 



Sm ' 6m 



nh),i~2L-T~'Qik2) 



For the quantities on the discrete spectrum one can notice first that 



Then 



Sm 



-0, + (_2AO-^ 
dm dm 

= (-2A;)^-"(n - l)Xict>iPi + {^2Xiy-"Qi. 



{Pl,{n),4>p,{n)}{n) 



Sm ' Sm 



% (-2L-i)"-i(-2Ap)i-" ((n - l)Xpcf>pPp + Qp) 



{n - l)Ap0p 



Pi,Pp 



(205) 



{P;,(n)'Pp,(n)}(n) 



— (-2L"i)"-i^ 
Sm ' Sm 



{<t>lXn)i<t>p,{n)}{n) 



= (-2A0^ 



(n - l)Xi(f>iPi + Qi, (n - l)Xp4>pPp + Qp 



(-2AO'""((n-l)Az0/ Pi,Qp +in-l)Xpcj)p 
(-2A0'""(n - l)(Aj0, - Xp(j)p)5ip = 0. 



□ 



10 Inverse scattering transform 

Inverse scattering method for the hierarchy (I17ip is the same as the one for the 
CH equation [^Oj. The only difference is the time-dependence of the scattering 
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data (|166p - (|168p . For example, the inverse scattering is simphfied in the 
important case of the so-called reflectionless potentials, when the scattering 
data is confined to the case Ti^{k) = for all real k. This class of potentials 
corresponds to the iV-soliton solutions of the CH hierarchy. In this case b{k) = 
and |a(fc)| — 1 and icLp is real: 



^ 1 + 2Kr. 



^e-^^li^^ -here « = E 1^ (1 



T^r) Kr) ~\~ K' 



Thus, ihp has the same sign as 6„, and therefore = bn/(iap) > 0. The time 
evolution of C+ is dTTO)) 

C+it) = C7+(0) exp ( - n,,niX-')t). (206) 

The iV-soliton solution is [50] 

uix,t) = '^ l\xp(^-\x-g{^, t)\)r^g^'{t m-co, (207) 
where g{^,t) can be expressed through the scattering data as 

g{U) ^ In / (l-E +1/2 (208) 

^ n,p ^ ' 

with 



^pn[^,t] = Spn + ■ 

Kp -\- Kyi 

For the peakon solutions {uj — 0) the dependence on the scattering data is 
also known E] . 

The CH multi-soliton solutions also appear in several works [HI [THl EH [Ml 
[eg [IQI [ni [111 [66]. The Darboux transform for the CH equation is obtained 
in [71] . The construction of multi-soliton and multi-positon solutions using the 
Darboux/Backlund transform is presented in [i^ I50j . 



11 Conclusions 

In this paper the Inverse Scattering Transform for the CH hierarchy is inter- 
preted as a Generalized Fourier Transform. The generalized exponents are the 
squares of the eigenf unctions of the associated spectral problem. Apparently 
the CH hierarchy is well defined only if q{x,0) = m{x,0) -I- w > 0. The only 
exception is the CH equation itself. The situation for CH when the condition 
q{x, 0) > on the initial data does not hold is more complicated and requires 
separate analysis |S21 [3 [H [13 [IH] • Throughout this work the solutions u{x,t) 
are confined to be functions in the Schwartz class, uj > 0. The inverse scattering 
is outlined in detail in [20] . 

The spectral problem ^ is gauge equivalent to a standard Sturm-Liouville 
problem, well known from the KdV hierarchy 

-(j)yy + U{y)(t) = /i=--^--A, 
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Uiv) - ^ + 9|«M_M^_± (210) 

Note that (|209p leads to two possible expressions for the change of the variables 
in the Liouville transformation: 



y — \/ujx + / {y/ q{x') — y/uj)dx' + const, (211) 



y = y/ujx+ i {y/qix') - y/u;)dx' + const. (212) 
These two possibilities, ()21ip . ()212p are only consistent iff 



(■\/ (7(x) — '\/aJ)dx = constant, 

which is always the case, since the integral under question is (up to a multiplier) 
the Casimir function a ((43|) : see some details in j61j . 

The matching of the CH hierarchy to KdV hierarchy requires solving the 
Ermakov-Pinney equation (|210p [HI US], which is not straightforward. One can 
eventually obtain a solution in parametric form [S31 [S3] , see also [70 , 71, 72 , 66] . 
The analytic properties of the eigenfunctions and especially their asymptotics 
for fc ^ oo in these two cases are substantially different [20], e.g. compare (j40)) . 
(|4ip with the well known results for the standard Stourm Liouville problem 
gtkx^^ + ...). This influences also the dispersion relation for the transmission 
coefficient. Thus the matching of the 1ST for these two cases is not automatic. 
This alternative approach relies on several implicit equations and is considerably 
less transparent than the approach adopted here. 

We have also excluded the possibility of 'creation' or 'death' of solitons, i.e. 
an appearance of a new discrete eigenvalue as a result of an infinitesimal change 
Sm in m. For the KdV equation this problem is addressed e.g. in [5^ . 

The behavior of the scattering data at A: = is also an important question. 
In our analysis we implicitly assumed that the Wronskian 

W{f+{x,k),r{x,k))U=o^O. 

Then a{k) at fc = has a singularity of type cf. (fTT)) . However it is possible 
that W{f^{x, k), f~{x, k))\x=n — 0, and then a{k) is not singular at fc = 0. To 
investigate the behavior of the quantities at fc = one can proceed as in [36j . 
There is a basis in the space of eigenfunctions of the spectral problem, which 
can be chosen as fi{x) = f~^{x,0) and f2{x) — —if^{x,0). The asymptotics 
are 



fi{x) — > 1, X oo; fi{x) Ax + B, x — oo; 
f2{x) — > x, a; — > oo; f2{x) —f Cx + D, x ^ — oo, 

where A, B, C, D do not depend on x, e.g. 
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1 f°° 

A = — m{x)fi{x)dx, 
1 f°° 

B = 1 / xm(x') fi (a;)da; 



1 /"^ 

C = 1 + — / rn{x)f2ix)dx. (213) 

The Wronskian has the same value at — cxd and +cxd, therefore 

BC-AD^l. (214) 
The behavior of the scattering data is 

a{k) = ^ + ^^ + "(1)' (215) 

b{k) = -± + ^—l + o{l). (216) 

Note that A is an integral of motion (cf. [HIllHI) as well as B + C since a{k) 
does not depend on i. If A = , the singularity disappears. Then, from (|214p 
it follows that EC = 1 (and B and C are integrals of motion by themselves in 
this case - since both B + C and BC are), i.e. if C = B~^ — ae^ , a = ±1, then 

a(0) = a cosh /3, 6(0) a sinh /3, 7^(0) = tanh /3. 

This situation, although exceptional is the one in which the purely solitonic 
case is allowed: when ^ = 0, 7^(0) = (if A 7^ 0, 7^(0) = -1, see (I^V 

The Poisson bracket ^ is defined through variations with respect to m. We 
established that among these variations there are some that do not vanish in the 
asymptotic limit — > cxd, such as for example the variations of the scattering 
data. This fact (which is not related to the smoothness or rate of decay of m{x)) 
is related to the presence (in general) of poles of a{k) and h{k) dX k — 0. More 
careful analysis [U [21 [37] leads to a modification in the definition of the Poisson 
bracket by additional terms, when the behavior of variations like 6a{k)/6m, 
db{k)/dm is considered at = 0. 
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